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Abstract 

We consider axion-mediated dark matter models motivated by 
Fermi gamma ray line at 130 GeV, where anomaly interactions 
of an axion-like scalar mediate a singlet Dirac fermion dark mat- 
ter (DM) to electroweak gauge bosons. In these models, extra 
vector-like leptons generate anomaly interactions for the axion 
and can also couple to the SM Higgs boson to modify the Higgs- 
to-diphoton rate. We can distinguish models by the branching 
fraction of the DM annihilation into a photon pair, favoring the 
model with a triplet fermion. From the condition that the lighter 
charged extra lepton must be heavier than dark matter for no tree- 
level DM annihilations, we also show that the ratio of Higgs-to- 
diphoton rate to the SM value is constrained by vacuum stability 
to 1.4(1.5) for the cutoff scale of 10(1) TeV. 



1 Introduction 



Both ATLAS and CMS [1] have recently discovered a Higgs-like boson with mass 125 — 
126 GeV with high significances of 5.9a and 5.0<r, respectively. Although the overall signal 
strength is consistent with what we expected from the SM Higgs boson, the Higgs-to- 
diphoton channel shows a larger signal strength than the Standard Model (SM) value. 
Since the Higgs decay rate into a photon pair is induced by loops in the SM, there might 
be a room for new particles to play a role in modifying the Higgs couplings to a photon 
pair without having been ruled out by collider limits and electroweak precision test yet. 

On the other hand, from the Fermi LAT data [2H5], there is an interesting observation 
of the gamma-ray line at 130 GeV [6], which might be interpreted as the signal of dark 
matter annihilating into monochromatic photon(s) at the galactic center [6H9]. There 
are a lot of recent activities in building the dark matter models [T0T[T2"] for explaining the 
Fermi gamma-ray line^with satisfying the bounds on other annihilation channels jUEHJdS]. 
Seemingly unrelated two phenomena with di-photons from the Higgs boson decay and those 
from dark matter annihilation could find the same origin beyond the SM. This is the main 
topic of this paper, addressing the interplay between dark matter and Higgs signals in 
constraining the model parameters and the signal strength expected from the models. 

We consider a microscopic theory of axion- mediated dark matter [TT], which has been 
recently proposed by ourselves to accommodate the Fermi gamma-ray line with a Dirac 
fermion dark matter. To that purpose, we introduce extra vector-like leptons to gener- 
ate the Peccei-Quinn(PQ) anomaly interactions of an axion-like pseudo-scalar to the SM 
electroweak gauge bosons. Then, the axion can mediate dark matter to annihilate into 
monochromatic photon(s) with the branching fraction depending on the representations of 
the extra leptons. Since we need a large coupling of the extra leptons to the Higgs boson 
to enhance the Higgs-to-diphoton rate [TBTfl9] . the minimal content of extra leptons that 
we consider contains: (1) one vector-like pairs of doublets and a vector-like pair of singlets 
or (2) a vector-like pair of doublets and a triplet. We need to introduce two Higgs doublets 
to allow for the Higgs Yukawa couplings to the extra leptons. However, we focus on the 
decoupling limit of extra Higgs bosons where only the SM Higgs boson is kept in the low 
energy theory. 

We show that Fermi gamma-ray line can be explained by the DM annihilation into 
a photon pair in both models. But, the two models are distinguishable by the different 
branching fractions of the DM annihilation cross section into a photon pair. When the 
lighter extra charged lepton state is smaller than around 200 GeV, the branching fraction 
of two-photon line increases significantly, as compared to the value determined by anomalies 
in the decoupling limit of extra lepton masses. In the singlet models, the branching fraction 
of two-photon line is large so extra annihilation channels with the scalar partner of the 
axion are needed. On the other hand, in the triplet model, the branching fraction of two- 

1 See also other interpretations in Ref. [13]. 
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photon line is consistent with the correct Fermi gamma-ray line and thermal relic density. 

We also show that the two models with extra leptons lead to the same Higgs-to-diphoton 
rate for the same charged fermion masses and couplings because the extra charged leptons 
are the same. However, vacuum stability conditions become stronger for the triplet model 
in the region with the enhanced Higgs-to-photon rate. In order to make the loop-induced 
DM annihilation into monochromatic photon(s) sizable for the Fermi gamma-ray line, we 
need to suppress or forbid the tree-level DM annihilation into a pair of the lighter state 
of extra charged leptons. Due to the nature of the axion coupling to the extra leptons, 
the phase-space suppression of the tree-level annihilation channel is not significant near 
where the lighter charged lepton is lighter than dark matter. Therefore, we need to forbid 
the tree-level channel kinematically by taking the lighter charged lepton mass to be larger 
than dark matter mass, 130 GeV. In turn, combined with vacuum stability conditions, we 
obtain the bounds on the ratio of Higgs-to-diphoton rate to the SM value as i? 77 < 1.4(1.5) 
for the cutoff scale of 10(1) TeV. 

The paper is organized as follows. We first present models with extra leptons and 
their interactions to the axion and the Higgs boson for the extended Higgs sector with 
two Higgs doublets and one complex singlet scalar. Then, we discuss the implications of 
extra leptons on the DM annihilation cross section into monochromatic photon(s) and the 
Higgs-to-diphoton rate, impose the vacuum stability and perturbativity conditions, and 
electroweak precision constraints, and emphasize the interplay between dark matter and 
Higgs signals. Finally, the conclusions are drawn. 

2 Models with extra leptons 

We consider the extension of the SM with renormalizable interactions for axion-like scalar 
field a, an extra vector-like lepton / and a singlet Dirac fermion dark matter x, as follows, 



where _D M is the covariant derivative with respect to the SM electroweak gauge group. 

When the extra leptons carry electromagnetic charges, they induce the effective inter- 
action between the axion and two photons by triangle one-loops. Consequently, due to the 
direct axion coupling to the fermion DM, monochromatic photons can be generated from 
the DM annihilations mediated by the axion. For the large extra lepton mass rrif — \/v s 
with v s being a singlet vev, PQ anomalies lead to the axion couplings to the electroweak 
gauge bosons, 
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m x xx ~ rrifff - -m 2 a a 2 + —i\ x a X ^X + i^afl 5 f 



(1) 




(2) 
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where c±,C2 are constant parameters fixed as c\ = Tr(qpQY 2 ) and c% = Tr(gpQ/(r)) with 
gpQ(y) being PQ charge (hypercharge) and Z(r) being the Dynkin index of representation 
r of the extra lepton under the SU(2)l- For a pair of extra lepton doublets like Higgsinos 
in the MSSM, we get c\ = C2 = 1 for PQ charges normalized to 1. 

Furthemore, the extra vector-like lepton can have the effective interaction with the SM 
Higgs H after heavy fermions are integrated out, 

Cn,es = jj^Hff. (3) 

In the presence of multiple vector-like leptons fi(i = 1, 2, ■ • • , N), the Higgs to di-photon 
decay rate can be enhanced or reduced, depending on the interference between the SM 
contribution and the new contributions with couplings Cf v 

2.1 Model I: vector-like doublet and singlet with PQ charges 

For concreteness, we consider PQ symmetry for a Dirac fermion dark matter that couples 
to the SM through the extra vector-like leptons. We introduce a Dirac singlet dark matter 
fermion composed of x, X\ an d extra vector-like lepton doublets, U, with I4 = (z/4, e^) T and 
h = (£4, va) t 1 and extra vector-like lepton singlets, e^, e%. A complex scalar, S, mediates 
between dark matter and extra vector-like leptons while two Higgs doublets, H d , H u , are 
required to obtain the Higgs Yukawa couplings to extra leptons being consistent with PQ 
symmetry. 

The interaction terms for dark matter and extra leptons are written in terms of two- 
component Weyl spinors as follows, 

- Yukawa = ^ X S XX + hShk + KSe\e\ + ViB d l^e\ - yiHJ A e c 4 + h.c. (4) 

Here we note that the vector-like lepton doublet and singlet have the same PQ charges and 
the Higgs doublets have the same PQ charges as the one of the singlet S. In the presence 
of PQ symmetry, the Yukawa couplings for quarks and charged leptons in the SM can be 
also written by assigning the PQ charges appropriately. 

For dimension-5 Majorana neutrino mass terms, jj(liH u )(ljH u ), where M is the cutoff, 
we need to choose nonzero PQ charges for lepton doublets but zero PQ charges for lepton 
singlets. For Majorana neutrino case, the mixing term between the SM charged leptons 
and the extra charged leptons such as HJ^, H d lie c A) are forbidden by PQ symmetry. 
Furthermore, the mixing dimension-5 terms between the extra heavy neutrino and the SM 
neutrinos are not allowed by PQ symmetry, so the lighter charged lepton would be stable 
in the minimal case. The PQ charges and Z2 parities in Majorana neutrino case are given 
in Table 1. 

On the other hand, if the PQ charges of lepton singlets are nonzero, only Dirac-type 
neutrino masses are allowed with right-handed neutrinos iV?. For Dirac neutrino case, 
the mixing term between the SM leptons and the extra vector-like leptons are allowed 
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Table 1: PQ charges and Z 2 parities for Majorana neutrino case 



by the renormalizable Yukawa couplings so there should be an additional symmetry to 
protect large flavor violations in the lepton sector. The PQ charges and Z 2 parities in 
Dirac neutrino case are given in Table 2. 
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Table 2: PQ charges and Z2 parities for Dirac neutrino case 



We note that in both Majorana and Dirac cases, Yukawa couplings are invariant under 
a shift in PQ charges, which reflects another accidental ^(1)^ global symmetry. 

2.1.1 Axion and Higgs couplings to extra leptons 

The PQ invariant potential for two Higgs doublets only is given by 

Vo{H u , H d ) = ii\\H d \ 2 + fil\H u \ 2 + -XilH^ + -X 2 \H U \ 4: + X 3 \H u \ 2 \H d \ 2 + X 4 \H u H d \ 2 . (5) 

The above Higgs potential V possesses an additional global U{1)h symmetry. Then, 
keeping the U{1)h in the singlet-Higgs couplings, we get the PQ and U(1)h invariant 
scalar potential for the singlet and the Higgs doublets, 

V(S, H u , H d ) = A S |S| 4 + m 2 \S\ 2 + A^H^I 2 + X HdS \S\ 2 \H d \ 2 + V {H U , H d ). (6) 

Thus, after PQ and U(1)h symmetries are broken spontaneously, there are two massless 
axions. If the U(l) symmetries are broken by the following soft mass terms, 

AV = ^m ,2 S 2 - ti 2 3 H u H d + h.c, (7) 

the massless axions obtain masses but there is no mixing between the S and Higgs axions. 
In both Majorana and Dirac neutrinos, the SM quarks are also charged under the PQ sym- 
metry, so they could contribute to QCD anomalies through the Higgs Yukawa couplings. 
However, when two massless axions are lifted by the above U(l) breaking soft masses, 
there is no mixing between the Higgs pseudo-scalar and the singlet axion. So, the axion 
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mediates dark matter interactions to electroweak gauge bosons through the couplings to 
extra vector-like leptons but it does not couple to the SM fermions or gluons directly. 

In the limit of decoupling the heavy additional Higgses, we can get the following effective 
potential containing the lighter CP-even Higgs (in unitary gauge) and the singlet only, 

V eS = \ S \S\ 4 + /4|S| 2 + (im'V + h.c.) + \x hS \S\ 2 h 2 + \\ h h* + \rn\h 2 (8) 

where Xhs, A/i,m 2 are the effective Higgs parameters written in terms of the input param- 
eters in V tot . For simplicity, henceforth we consider the effective theory with the heavy 
Higgses being integrated out. In this case, if there is no mixing between the lighter Higgs 
and the singlet, the lighter Higgs can be just like the SM Higgs. 

After the singlet scalar and the Higgs doublet develop VEVs as (S) = v s , Hd = 
■^{vd, 0) T and H u = ^(0, v u ) T with v 2 = v\ + v\, we get the mass terms, 

- £ mass = m x xx + mi(u A U4 + e 4 e 4 ) + m e e\e\ + m x e 4 el + m x e A e\ + h.c. (9) 

with m x = ^X x v s , m = Mv s , m e = X e v s , m x = -^yiv d and m x = ^ViV u . Then, the 
charged fermion mass matrix is 

(W 1"' do) 



v m x m e J \e\ / 
Thus, the mass matrix squared is given by 

M\M f =( m " + < m f* + 2 m ° m * ) . (ii) 
7 y rnim x + m e m x m e + m x J 

Then, after diagonalizing the squared matrix by a unitary transform with V, 

V ] M ] MV = M 2 D (12) 

with 



(13) 



y = f COs6 f sin #/ \ 

\ — sin 8f cosOf J 
we get the mass eigenvalues and the mixing angle are 

< 2 = l(r^ + ml + ml + mlTy/ K + < - ml - m')* + ±{m x m x + m e m x y) , (14) 

= _ 2(m l fh+mm x ) _ 

mi — mi 

h J2 

Accordingly, the Weyl fermion pair (e4,e1) T transforms as 

(SMa)- (16) 
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By performing a bi-unitary transform WAi V = Mp with 

(^-{mi cosOf — m x sinOf) ^-(m; sin 6^ + m x cosOf) 
^-j-(m x cosOf — m e sin #j) -^-{m x sin Of + m e cos Of) 

we can diagonalize the charged fermion mass terms in eq. ({TO]) as follows, 

m h fxfi+m h f 2 h 
Then, another pair of Weyl fermions, (§4, el) T , transforms as 



(17) 



U* 



h 
h 



(19) 



Now we expand the singlets and Higgs doublets around the VEVs by S = (v s + s + 
ia)/y/2, Hd = ^(^d + h d ,0) T and H u = -^(0,v u + h u ) T . In the decoupling limit for the 
extra Higgs fields, the lightest CP-even Higgs becomes the SM Higgs while the lighter axion 
coming from S remains in the effective theory. Then, from the results in appendix A, we 
get the Yukawa couplings for the axion a and the Higgs h, 



,m y _ 5 .mi _ 5 
-1 — - ax'J X ~ 1 — av l v 
v s v s 

-ia (\ 1 F 1 -f 5 F 1 + A 2 F 27 5 F 2 



(A 3 - A 4 )(F 2 i 7 \ - F 1 F 2 ) - (A 3 + A 4 )(F 2T 5 J P 1 + F l7 5 F 2 ) 



1 

+ 2 m _ 

-yihFxFi - y 2 hF 2 F 2 

1 ' f(y 3 + Va) {F 2 F t + F,F 2 ) + (-y 3 + y A )(F 2l 5 F 1 - F ll5 F 2 ) 



-h 



(20) 



where x = (x> X^)i v = (^4; v \)-> Fi = fi) T an d -^2 = {h, f 2 ) T are the mass eigenstates 
in Dirac spinors. For simplicity, we take mi = m e , m x = m x for which Of 



\ and 



m /i 2 = m i T m x- Then, the Yukawa couplings are 



Ai 
Vi 



V s 

m x 
v 



In this case, eq. ( 1201) becomes 

— £int 



A2, A3 — A4 — 0, 
: -V2, 2/3 = 2/4 = 0. 

.m 



(21) 
(22) 



-1 — a xi X ~ 1 — a 1/7 v 

Vs Vs 



-zAia {F ll b F l + F 27 5 F 2 ) - y x h {F X F X - F 2 F 2 ). 



(23) 
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On the other hand, using the results in appendix B, we get the electroweak interactions of 
the extra leptons, 



''-gauge 2 V 



^(Fi + F 2 ) + h.c. 

v-fv + v e (F 1 ^F 1 + F 2 ^F 2 ) + a e (F 1 ^F 2 + F?rfF x ) 
-eA fl (F 1 YF 1 + F 2 YF 2 ) (24) 



2 sin 9 W cos 6w M 



where v e = —\ + 2sin 2 6*vy and a e = — \. We note that there are also Yukawa couplings 
between the CP-even singlet scalar and dark matter/extra leptons. If they mix with the 
Higgs, they can affect the dark matter annihilation as well as the effective Higgs couplings 
to the SM. Here, we assume that the mixing between the CP-even singlet scalar and the 
Higgs is small enough not to affect the Higgs production cross section at the LHC. 

2.1.2 Dirac neutrino and SU(3)w unification 

For Dirac neutrino case, each pair of the SM lepton doublet and singlet also carries the 
same PQ charges as for the extra leptons, so the case is consistent with the SU(3)w 
unification (20] in which the electroweak gauge group SU(2) L x U(l)y is unified into a 
single one. For instance, in 5D SU (3)w electroweak unification model on S 1 /(Z 2 x Z' 2 ) [2Tj . 
SU(3)w is broken to SU{2) L x U(l)y due to the orbifold boundary conditions, resulting 
in two fixed points on the orbifold: y = where SU(3)w is unbroken and y = ^ with R 
the radius of the extra dimension where SU(3)w is broken to SU{2) L x U(l)y. Thus, the 
SM leptons and extra leptons, that are a full SU(3)w triplet, are located at y = while 
the SM quarks, that are split multiplets of the SU(3)w, are located at y — For the 
realistic Yukawa couplings for the charged leptons with SU(3)w, the Higgs doublet must 
be extended to a bulk sextet, ifg = (Hd, Ht, Hs), including a doublet Hd with Y = —1/2, 
a triplet Ht with Y = — 1 and a singlet Hs with Y = —2, instead of an SU(3)w Higgs 
triplet. There is one more bulk Higgs sextet, Hq = (H u , Ht, H$), which contains another 
Higgs doublet, H u . Then, the Higgs doublets survive as zero modes, giving rise to the 
Yukawa couplings for the SM quarks, leptons and extra leptons, while the extra Higgs 
states are projected out by orbifold boundary conditions. 

In 5D SU(3)w unification on the orbifold, the KK towers of Higgs sextets living in 
bulk contribute to the running of the electroweak gauge couplings, Under the assumption 
that the unified electroweak gauge coupling becomes strong, we get the prediction on the 
Weinberg angle at Mz with KK corrections 
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8tt 



M, ... M, 



sin 2 e w (M z ) = 0.25 - — a em B In — *- + B In — 1 (25) 



Mr. M Z i 



where a(M z ) = f- ^ 1/128, B = b g — b g >/3 and B = b g — b g/ /3. Here b's denote the 
beta function coefficients for the SM zero modes and b's denote those for the KK modes. 
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In the case with one pair of Higgs sextet and anti-sextet, the beta functions below the 
compactification scale are b g = 3,b g > = — y and B — Above the compactification 
scale, we get b g = — ~, b g > = — ^ and B = ||. Then, from the observed value with 
sin 2 9 w(Mz) ~ 0.231, we can determine the unification and compactification scales under 
the strong coupling assumption with M*/M c ~ 8tt 3 (100) as M* = 140(77) TeV and 1/R = 
1.7(2.4) TeV. 

In the case of the SU(3)w unification, PQ anomalies respect the 577(3) iy symmetry as 
follows, 

C a , sui3)w = C £^aF^F^ (26) 

2 

with a = with g being the 577(3) w gauge coupling and c = 1. Consequently, from 
the embedding of U{1)y into SU(3)w with Y = diag(— ~, — ~, 1) and = <7iF, we get 
the relation between the S77(3)w gauge coupling and the electroweak gauge couplings as 
go = \/3gi = 02 at tree level. Thus, we obtain the tree-level value, sin 2 9 W = /- 1 a = 0.25, 
at the 5?7(3)vk unification scale. Including the running of the gauge couplings below the 
unification scale, the electroweak anomaly couplings are given by the same formula (12"B"|) 
but in terms of electroweak gauge couplings, 

^-E^ a V p ( 27 ) 

i=l,2 s 

with C\ = 3c 2 . For an arbitrary vector-like pair of the 577(3) representations, the ratio 
is fixed to C\/c2 = 3. For a pair of SU(3)w triplets that we consider for minimality, we 
get c\ = 3 and c 2 = 1. The result is in contrast with the case where there are only a pair 
of extra lepton doublets like Higgsinos in MSSM for which c\ = c 2 = 1. Furthermore, we 
also note that the S77(3)w requires a simpler mass matrix of extra charged leptons by the 
relation, mi = m e . 



2.2 Model II: vector-like doublet and PQ-neutral singlet 

We assume that a vector-like doublet lepton has nonzero PQ charges but a vector-like 
singlet lepton has zero PQ charges. The corresponding interaction terms for dark matter 
and extra leptons are then the following, 

- C = \ x Sxx + \Skk + m e ele c 4 + yiH d he% - yiHJ A el + h.c. (28) 

where the vector-like singlet lepton has nonzero tree-level mass. 

For Majorana neutrino case, the SM lepton doublets have the same PQ charges as the 
extra lepton doublets, so the mixing term between the SM leptons and the extra vector- 
like leptons are allowed by the renormalizable Yukawa couplings. The PQ charges and Z 2 
parities are given in Table 3. 
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Table 3: PQ charges and Z 2 parities for Majorana neutrino case 



For Dirac neutrino case, we can forbid the Majorana neutrino mass term by choosing 
the PQ charges of the SM lepton doublets to be different from the extra lepton doublets. 
In this case, there is no mixing between the extra leptons and the SM leptons. The PQ 
charges and Z 2 parities including the right-handed neutrinos N£ are given in Table 4. 
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Table 4: PQ charges and Z 2 parities for Dirac neutrino case 



In the model with PQ-neutral extra singlets, the electroweak anomalies come only from 
the vector-like doublet lepton with anomaly coefficients, C\ — c% — 1. Although the Higgs 
to diphoton decay rate is not changed as compared to the case with extra singlets with 
nonzero PQ charges, we get the branching fraction of DM annihilation into a photon pair 
to be about 14% in the limit of large extra lepton masses as will be discussed in later 
section. 

The Higgs interactions to extra leptons are the same as in Model I. But, the axion 
interactions are different from those in Model I because the singlet leptons do not couple 
to the axion. For mi = m e and m x = m x , the axion interactions in eq. (T20]) become 

Ai = A 2 = A 3 = A 4 = p-. (29) 

2.3 Model III: vector-like doublet and triplet with PQ charges 

When the extra vector-like singlet fermion is replaced by a triplet fermion T with zero 
hypercharge, the Higgs Yukawa couplings in eq. (jl]) become 

- Yukawa = h ViH d kT - ViHJ^T + h.c. (30) 

Then, as in the model with vector-like doublet and singlet leptons, there are two charged 
Dirac fermions with masses m/^m^. The effective Yukawa couplings for the axion and 
the lighter CP-even Higgs are obtained similarly as in eq. (!20l . We can assign the same 
PQ charge and Z 2 parity for the triplet fermion as for the vector-like singlet fermion. 
Consequently, as far as the Higgs to diphoton rate is concerned, the triplet model is similar 
to Model I. 
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The triplet model has crucial differences from Model I. First, vacuum stability bounds 
are stronger in the triplet model due to a more negative contribution of the triplet fermion 
in the RG equation of the Higgs quartic coupling [17] . Second, in the triplet model, the 
coefficients of the axion interactions to electroweak gauge bosons are given by c% — 1 and 
C2 = 3, in the limit of large extra lepton masses, in contrast to c\ = 3 and C2 = 1 in Model 
I. Therefore, the intensity of the two-photon line is comparable to the one-photon line. In 
this case, the Fermi gamma-ray line can be still explained by the two-photon line, when 
dark matter mass is 130 GeV [12] . 



3 Interplay between dark matter and Higgs signal 

In this section, we consider various constraints on the model with vector-like leptons, from 
the DM annihilation cross section into monochromatic photons, the Higgs to diphoton 
decay rate and electroweak precision data as well as cosmology and collider bounds. 

3.1 DM annihilation into monochromatic photon(s) 

From the interaction between the axion and the vector-like charged extra lepton in eq. (CQ), 
we get the analytic expression for the amplitude for the axion decay into a pair of photons, 
M = e*^{k l )e* 1/ {k 2 )M^, with 

■Mau = i — ~^vptiu^2^\ (31) 

7r rrif 

where 77 = Am 2 t / m 2 a and 

A x {t) = rarcsin 2 (l/v^). (32) 
In the limit of r ^> 1, A\{r) « 1, we obtain the following, 

■M-uv ~ i ■ tvpiMjk^ki . (33) 
7T m f 

Then, we can approximate the axion coupling to a pair of photons as 

•^077 Cd^ttF^uF^ 1 (34) 

with 

A/OJem 1 



, (35) 

For the axion effective interactions to electroweak gauge bosons, we can take the sim- 
ilar limits of large fermion masses running in loops or equivalently compute the anomaly 
coefficients. Then, using the axion effective interactions, from Ref . [11] , we obtain the DM 
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Model I 


Model II 


Model III 


(ci,c 2 ) 


(3,1) 


(1,1) 


(1,3) 


Br(xx -> 77) 


40% 


14% 


6.5% 


Br(xx -> WW) 


44% 


62% 


65% 


Br(xx -> ZZ) 


16% 


16% 


15% 


(r,R) 


(1.15 x 10- 3 ,0.56) 


(0.27, 1.77) 


(1.01,2.51) 



Table 5: Branching fractions for DM annihilation cross sections in the decoupling limit of 
extra leptons. Here, r = (av) Zj /(2(av) 11 ) and R = {av) ww /{2(av) 11 + (av) Zj ). 



annihilation cross sections into 77, Zj, WW and ZZ, respectively, in terms of the anomaly 
coefficients, ci,C2, in eq. 

|A I 2 

(cru) 77 = x (ciai cos 2 6*vy + c 2 a 2 sin 2 9 W ) 2 S, (36) 

(av) Zl = 47^^(c2«2- Cl a0 2 sin 2 (2M5(l-^) 3 , (37) 

{av)ww ~ 102W S S V-^J ' (38) 

I \ 12 , v 3/2 

= ^T^TT (c 2 a 2 cos 2 e w + c l a 1 sm 2 9 w ) 2 S(l f) ' (39) 



where 



s 2 



S = 7 2 r 2 , r2 2 : * - 4m 2 . (40) 

(s - mi) 2 + Y 2 a m 2 a 



Then, the ratio of the intensities of two-photon to one-photon lines is 

r _ {vyhj (c 2 a 2 - ciai) 2 sin 2 (2^) / M|\3 

~~ 2((7v)~~ 4fci«i cos 2 Ow + coao sin 2 9 W ) 2 v 4m 2 / 



x 



Therefore, we get the following results with m x = 130 GeV: in Model I with c\ — 3, c 2 — 1, 
r = 1.15 x 10~ 3 ; in Model II with Ci = c 2 = 1, r = 0.27; in Model III with d = 1, c 2 = 3, 
r = 1.01. In Model III, even if the intensity of two-photon line is comparable to the one of 
one-photon line, the best fit with two lines occurs for m x = 130 GeV as in the two former 
models [12] . In Table we summarize the branching fractions for DM annihilation cross 
sections in the decoupling limit of extra leptons. The fraction of the continuum photons, 
R = {av) ww /{2{av) 11 + (av) Z7 ), is also shown in the same table and we find that the 
continuum photons produced from the W^ll^-channel would be consistent with the line 
spectrum in all the three models |15j . 

Now we consider the effect of the finite extra lepton masses in the DM annihilation 
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cross sections. The amplitude for the DM annihilation into a pair of photons is 



M 



xx-m 



(M 



s - ml- iY a m a 



{M a - 



■11) 



(42) 



where s is the center of mass energy squared and T a is the total decay width of the axion-like 
scalar. Then, the annihilation cross section is given by 



(av) 



1 



77 



167TS 



\M 



xx-m I 



(43) 



Here, the squared amplitude for a — > 77 is obtained from the result of the decay width of 
the axion to two photons at one loop, 



1 



a— >-yy 



32iim a 

\x f \W 



\M 



a— >77 1 



64vr 3 



" m ml 



Mrf) 



(44) 



Since the axion-like scalar can be off-shell in the s-channel, we need to replace m a in eq. (T44 
with s for the correct |A^ a -s>77| 2 in eq. (j4"3l . Thus we get 



\M a - 



•77 1 



l A /l ft cm 2 
27T 2 



Mr/) 



m f 



(45) 



where Tf is now replaced by Am'j/s. Therefore, plugging the above into eq. fj43|) with 
eq. (jl2l) . we obtain 



A X A / | 2 « 2 



(av) 



77 



512tt j 



fs - m 2 ) 2 + r 2 m 2 



Mrf) 



m f 



(46) 



In Model I, there are a pair of mass eigenstates of extra charged leptons with axion 
couplings, Ai, A2, which are given in (J2TJ) . So, we get the DM annihilation cross section to 
a photon pair as 



(av) 



IA l 2 a 2 

I^Xl "cm 



77 



5127T 3 (s — m 2 ) 2 + r 2 m 2 



Ai^i(n) Aiv4i(r 2 ) 



771/1 



m /2 



(47) 



with Tj = 4m 2 ./s. On the other hand, from the results of appendix C, the DM annihilation 
cross section to Z-boson and one photon is 



(av) 



Zj 



|A x | 2 a em Q; (-0.5 + 2sin 2 6" w ) 2 



1024tt s 



cos 2 6 W 



(s — m 2 ) 2 + r 2 m 2 



m| ^ 3 



-V 



4m 2 / 



x 



m /2 
12 



(48) 




m a (GeV) 



Figure 1: Parameter space of m a and A x where the annihilation cross section into a photon 
pair in Model I. Yellow region is between the central values of the cross section for Einasto 



and NFW dark matter profiles, 10 GeV c and 2 x 10 GeV c, respectively. We took 
Vi = Hi/ tan/3 for m x = fh x and rrih = 125 GeV for Higgs mass. We chose mf x — 140 GeV 
and m/ 2 = 460 GeV, which leads to _R 77 ~ 1.4 being consistent with vacuum stability at 
ITeV. 



where pf = 4m 2 . /m| and 



rp 



v 



T — p v 



jarcsin 2 (l/\/~~f^\ — arcsin 2 ( 1 jyfr) j . 



(49) 



Therefore, we get the ratio of two-photon to one-photon cross sections for Model I as 



(av) 



Z-y 



a(-0.5 + 2sin 2 6 



w , 



(av) 



ii 



2a em cos 2 8 W 



1 - 



m 2 7 \ 3 



AiAt(ri,pf )/m h + \ 2 A 4 (T 2 ,p2)/m h 



Ami, 

X, 



\ 1 A l (T 1 )/m h + \ 2 A l (T 2 )/m h 



(50) 

So, for sin 2 8yy(Mz) — 0.231, the overall factor of the ratio is very small so the one- 
photon channel is negligible as the factor containing the loop functions is of order one. 
Then, it is possible to obtain the gamma ray line at 130 GeV from the two-photon peak, 
when dark matter mass is 130 GeV. We note that the DM annihilation cross section into 
a pair of photons is required to be (av) 71 = (1.27 ± 0.32^28) x 10~ 27 cm 3 s _1 for the 
Einasto profile and (av)^ = (2.27 ± 0.57±^ 2 ) x lO^cmV 1 for the NFW profile, that is, 
Br(xx ~^ ll) — 4 — 8% for thermal dark matter [6]. Therefore, since the DM annihilation 
into a photon pair is loop-suppressed, we need to rely on a resonance effect near m a ~ 2m x 
to get the correct annihilation cross section [UJ. In Fig. [U we depict the parameter space 
of the axion mass m a versus dark matter coupling A x , satisfying the required annihilation 
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Cutoff 10TeV, tan(P)=1 
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Figure 2: Branching fractions of the DM annihilation cross sections into a photon pair, 
WW and ZZ, as a function of the lighter charged lepton mass in Model I. Gray region 
is for rrif 2 < lOTeV and red region is for rrif 2 < 2TeV. We varied over < A x < 1 and 
took yi = y\j tan/3 for m x = rh x and rrih = 125 GeV for Higgs mass. In red region, the 
annihilation cross section, (cri>) 77 = (1.27 — 2.27) x 10 _27 cm 3 s _1 , is obtained. Blue region 
is consistent with perturbativity, vacuum stability and EWPD within 2a, which will be 
discussed in later sections. 

cross section into a photon pair. We can see that the required tuning in the axion mass 
can be smaller for a larger dark matter coupling. For instance, for A x < 0.8, the axion 
mass can be deviated by \m a — 2m x | < 10 GeV from the resonance mass. 

Similarly, as given in appendix C, we also obtain the matrix elements for the other 
channels such as the intermediate axion going to ZZ and WW for Model I. The WW- 
channel is induced by loops containing the axion couplings to the heavier charged lepton 
and the extra neutrino. On the other hand, charged fermion leptons and extra neutrino 
couple to the axion with similar strength, in particular, with the same strength for the 
maximal mixing of charged leptons. Thus, for a small mass of the lighter charged fermion, 
the branching fraction of the IWT-channel is much smaller than the one of the channel 
with a photon pair, even if they are similar in the limit of large charged lepton masses. As 
shown in Fig. [21 constraints coming from perturbativity, vacuum stability at 10 TeV and 
EWPD within 2a are stronger than the one from Fermi gamma-ray line. In the parameter 
space consistent with all the constraints, the branching fraction of two-photon line depends 
on the lighter charged fermion mass: it becomes larger than 50% for < 200 GeV while 
it approaches 40% in the limit of 3> m x . Moreover, the WW channel is next-to- 
dominant and it becomes smaller than 40% for < 280 GeV and it approaches 44% in 
the large fermion mass limit. Now we recall the 77-channel should occupy about 4 — 8% 
branching fraction of thermal cross section for explaining the 130 GeV Fermi gamma-ray 
line, depending on the dark matter profile. Thus, in order to get the correct branching 
fraction of the 77-channel in Model I, there must be extra channels for DM annihilation 
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at freezeout with at least (av)x = (6.5 — 14)(av) a where (av) a is the annihilation cross 
section from axion mediation. The scalar partner of the singlet axion can give rise to 
extra DM annihilation channels into a pair of the SM particles through the Higgs-singlet 
coupling [TTJ [22]. These extra channels are p-wave suppressed so they can contribute 
to thermal cross section at freezeout while being temperature-suppressed at present [TTj . 
Furthermore, the Higgs-singlet coupling can be constrained by DM direct detection and 
Higgs production cross section at the LHC [TTj . 

On the other hand, in Model II, the vector-like doublet lepton couples to the axion 
but the vector-like singlet does not, so the 77 channel for DM annihilation has a smaller 
branching fraction than in Model I. For simplicity, in the maximal mixing case with mi = 
m e , m x = m x , from the axion interactions in eqs. (12"Tj) and (129]) . the DM annihilation 
cross section into a photon pair in Model II is reduced by a factor 4 as compared to the 
one in Model I. Furthermore, there are extra diagrams producing WW, Z7 from the DM 
annihilations due to the level-changing axion interactions, so the branching fraction of the 
77 channel in Model II is reduced roughly by a factor 3 to about 14% as compared to Model 
I, as expected from the anomalies. Therefore, the extra channels coming from the scalar 
partner of the axion should give a smaller contribution of at least (av)x = (0.8 — 2.7) (av) a 
at freezeout than in Model I. 

Finally, in Model III, instead of a vector-like singlet, a triplet is introduced with the 
Yukawa couplings to the axion and the Higgs. In this case, although the annihilation cross 
section into a photon pair is similar to half the one into Z7, the two-photon line can still 
explain the Fermi gamma-ray line at 130 GeV [12]. The DM annihilation cross section into 
a photon pair is about 6% in the limit of large extra lepton masses and it increases a bit 
when the lighter charged state is light. Therefore, in Model III, without extra annihilation 
channels at freezeout, the axion mediation only can accommodate Fermi gamma-ray line 
by dark matter with thermal cross section. 

3.2 Electroweak precision constraints 

With the improved value of the observed W-boson mass [23], the global electroweak pre- 
cision analysis with the Higgs mass at = 125 GeV gives rise to the constraint on the 
STU parameters as [23] 



S = 0.00 ±0.10, T = 0.02 ±0.11, U = 0.03 ± 0.09. 



(51) 



The contribution of the vector-like extra leptons to the T parameter 




is 



AT 



^[X)(l^| 2 + l C/ H| 2 )^+(^^)±2Re(FHC/ 1 * i )^(^|/0 



{\V 12 V U \ 2 + \U 12 U 11 \ 2 )9 + (y 1 ,y 2 ) - 2Re(V 12 V* 1 U* 12 U 11 )9_(y 1 ,y 2 ) 



(52) 
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Cutoff 1 0TeV, tan(P)=1 Cutoff 1 OTeV, tan<P)=3 




m f (GeV) m, (GeV) 



Figure 3: Parameter space of vs mj 2 for the singlet models with cutoff Mj = lOTeV. 
Solid lines show i? 77 = 1.0 — 1.7 with step 0.1 from bottom to top. The region above 
red solid(dashed) line is disfavored by EWPD at 68(95)% C.L. Blue region is excluded 
by perturbativity while gray region is excluded by vacuum metastability. tan/3 = 1(3) is 
chosen for the left(right) plot. We took yi = yij tan/3 for m x = m x and rrih = 125 GeV for 
Higgs mass. 



where yi = m 2 /m 2 Z) y\ = mf/m%, and 



n , i\ i 2ab a 

V+{a,b) = a + b r in 7; 

a — o b 

6-(a,b) = 2Vab[ —— In- -2). 

\a — b b 



(53) 
(54) 



In the case with mi = m e , m x = m x , the above contribution to the T parameter becomes 

1 



AT 



-^ + {yx,y 2 ) - -6-(yi,y 2 ) 



(55) 



As AT is more stringent than AS in the models |16] . it is enough to impose AT for the 
EWPD constraint in next section. 



3.3 Higgs to di-photon decay rate 

The heavy charged leptons, introduced to explain the Fermi gamma-ray line, can also lead 
to extra contributions to the decay of the SM-like Higgs into a photon pair through the 
Higgs couplings. In this section, we discuss the impact of EWPD constraints, perturbativity 



16 



Cutoff 1TeV, tan(P)=1 



500 

> 

2-400 



- 


V 


7*-'A/ : / 
/V/ "/ . 
































































I I I 


1 1 1 1 1 1 1 


Ml 


m f > m 

1 1 1 1 1 1 1 


U 

in 


1 1 r 



Cutoff 1TeV, tan(P)=3 



TOO 120 140 160 180 200 220 240 260 280 300 

m f (GeV) 



500 

> 
CD 

2,400 



- 
- 




/ : X** 




// / ^ 




















A-"' I !-- 






< 


>^ : 
















m j > rti f I 


j-rf 


I I I I 1 1 I 


1 1 1 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 



100 120 140 160 180 200 220 240 260 280 300 

m, (GeV) 



Figure 4: The same as in Fig. [3J but with cutoff Mj = 1 TeV. 



and vacuum stability conditions on the parameter space for which the Higgs diphoton signal 
is enhanced. Furthermore, we show how Fermi gamma-ray line restricts the parameter 
space further. 

We define the ratio of the Higgs production cross section times the branching fraction, 
yU 77 = (o-x^f r7 )g M • ^ ne re P or1ce d signal strengths for the Higgs-like events from ATLAS and 
CMS are the 7 following pQ, 



/i^ LAS = 1.90 ± 0.5, 



CMS 
M 77 



1.56 ±0.43, 



combi 
/i 77 



1.71 ±0.33. 



(56) 



The excess in the Higgs-to-diphoton signal might come from a statistical fluctuations or 
systematic uncertainties. Furthermore, it might well be that the QCD uncertainties were 
underestimated such that the deviations from the SM Higgs couplings could not be so sig- 
nificant [26]. Nonetheless, we investigate how much the Higgs diphoton signal is enhanced 
by the Higgs couplings to extra vector-like leptons in our models. 

In the presence of the effective Higgs coupling to charged fermions with Cf = —yf in 
eq. (jSJ), the ratio of the Higgs to di-photon decay rate with respect to the SM value becomes 



R 



77 



1 + 2// 



MTf) 



where N r 



m f A v {t w ) + N c Q 2 t A f (T t ) 

2 



(57) 



3 is the number of color, Q t = ±3 is the top quark electric charge in units of 
|e|, and Tj = Am 2 /ml, i = t,W, f. Below the WW threshold the loop functions are given 
by 



A v (x) 
A f {x) 



—x 
2x 2 



2x~ 2 + 3x- x + 3(2x" 1 - 



aT 1 + (a; -1 



(58) 
(59) 
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Cutoff 1 0TeV, tan(p)=1 Cutoff 1 0TeV, tan(P)=3 
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Figure 5: Parameter space of vs mj 2 for the triplet model with cutoff Mj = lOTeV. 
Solid lines show i? 77 = 1.0 — 1.7 with step 0.1 from bottom to top. Blue region is excluded 
by perturbativity while gray region is excluded by vacuum metastability. tan/3 = 1(3) is 
chosen for the left(right) plot. We took y\ — yi/ tan (3 for m x = rh x and = 125 GeV for 
Higgs mass. 



with f(x) = arcsin 2 -^ for x > 1. In the case with two charged Dirac fermions, the Higgs 
to di-photon decay rate is 



1 + v yh A f( T J/ m h + Vh A A T i)l m h 



A v {t w ) + N c Q 2 t A f (r t ) 



(60) 



We note that the above ratio can be compared to the experimental measure by i? 77 ~ /i 77 
when the singlet mixing to the SM Higgs boson is small enough. 

In Figs. [3] and HI in the case of the maximal mixing of extra charged fermions for 
Model I and II (the singlet models), we show the parameter space giving rise to the Higgs 
to diphoton decay rate larger than the SM value, depending on the cutoff scale Mj and 
tan/3 = v u /vd- We have overlaid EWPD constraints at 68% and 95% C.L., perturbativity 
and vacuum metastability by using the renormalization group equations in appendix D. 
We didn't include the coupling between the Higgs and the scalar partner of the axion in 
our analysis, although it could improve vacuum stability a bit if sizable. Furthermore, we 
assume that the axion couplings to dark matter and extra leptons are perturbative at least 
until the vacuum stability scale of the SM sector. As far as EWPD constraints at 95% C. 
L. are concerned, the Higgs to di-photon decay rate can be enhanced up to i? 77 ~ 1.6 for 
mf x < 115 GeV and mj 2 < 450 GeV. However, vacuum stability bound is the strongest 
constraint for the case with a large mass splitting so the Higgs to di-photon decay rate can 
be enhanced to at most i? 77 ~ 1.4 at < 110 GeV and mj 2 < 450 GeV for Mj = 10 TeV 
and i? 77 ~ 1.6 at < 110 GeV and rrif 2 < 360 GeV for Mi = 1 TeV. Consequently, in 
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Cutoff 1 TeV, tan(P)=1 Cutoff 1 TeV, tan(P)=3 
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Figure 6: The same as in Fig. |5l but with cutoff Mj = 1 TeV. 



order to get the Higgs to diphoton rate to be compatible with the LHC within la, we need 
a lighter charged extra lepton to be around 110 GeV. 

In Figs. |5] and El we show the parameter space for Model III (the triplet model) with 
perturbativity and vacuum metastability. Due to the more stringent vacuum stability 
bounds, we get at most i? 77 ~ 1.3 at < 120 GeV and rrif 2 < 360 GeV for Mj = 10 TeV 
and i? 77 ~ 1.5 at m fl < 115 GeV and m h < 440 GeV for M 7 = 1 TeV. 





Model I 


Model II 


Model III 


(ci,c 2 ) 


(3,1) 


(1,1) 


(1,3) 


Br(xx -> 77) 


>40% 


>14% 


>6% 




<l-5 


<l-5 


<l-4 



Table 6: Summary of the model predictions 



From the Fermi gamma-ray line, we require the tree-level DM annihilation into a pair 
of lighter charged leptons not to contribute, because otherwise the tree-level channel would 
easily dominate the DM annihilation into a photon pair. Then, we need the lighter charged 
extra lepton to be larger than dark matter mass, m x = 130 GeV, in both singlet and triplet 
models. Therefore, for > 130 GeV, in combination with vacuum stability conditions, 
the Higgs-to-diphoton rate in the singlet (triplet) model is bounded by i? 77 < 1.4(1.3) for 
the cutoff at Mj = 10 TeV and i? 77 < 1.5(1.4) for the cutoff at Mi = 1 TeV. The model 
predictions for the DM annihilation cross section into a pair and the Higgs-to-diphoton 
decay rate are summarized in Table [6j 
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3.4 Cosmology and collider bounds on extra leptons 

If extra charged leptons do not couple to the SM leptons directly as is the case in our 
minimal models, the extra lighter charged lepton would be stable. Then, it would be 
dangerous for the successful Big Bang nucleosynthesis(BBN) and the collider limits on 
stable charged particles. 

First, from the success of the BBN, we could put a bound on the lifetime of the lighter 
charged lepton to be smaller than about 5 x 10 3 sec, similarly to the case of stay NLSP 
in the MSSM [23], by considering the bound state effects of charged particle during BBN 
epoch. Furthermore, if the lighter charged lepton is stable, there are strong bounds on 
the mass of stable charged states at the LHC [29] so that one could not get a large Higgs 
to diphoton decay rate. Therefore, we need to extend the models p^6l[T7] . First, we can 
introduce a mixing of the extra doublet neutrino with an SM singlet neutrino. In this 
case, the lighter charged fermion can decay into VT-boson and the SM singlet neutrino 
sufficiently fast to satisfy the BBN and collider bounds. Second, extra charged lepton can 
mix with the SM charged leptons by the renormalizable Yukawa couplings such that the 
lighter charged lepton decays by f\ —> Zl and fi — > Wv. In this case, however, one has 
to take the sufficiently small mixing terms to satisfy the constraints coming from flavor 
violation such as fi — > e7 and r — > //7 and the Z-width. Requiring the mixing angle 
I Ua 4 1 < 1CT 2 (IT] , we can satisfy those constraints and make the lighter charged lepton 
decay prompt. 

4 Conclusions 

We have considered the possibility that extra vector-like leptons enhance the Higgs-to- 
diphoton rate and generate anomaly interactions in axion-mediated dark matter models. 
In order to get a large DM annihilation cross section to monochromatic photons as hinted 
from the Fermi gamma-ray line, the tree-level channels for DM annihilations must be 
suppressed or forbidden, constraining the dark matter models. In our models, the masses 
of the extra leptons that couple to the axion mediator must be larger than dark matter 
mass. Furthermore, dark matter with thermal cross section restricts the branching fraction 
of the DM annihilation cross section into a photon pair, favoring the triplet model in which 
the SU(2)l anomalies are larger than in the singlet models. 

The enhancement of the Higgs-to-diphoton rate requires new large couplings between 
the Higgs and new light charged states, so it can be comparable to the signal strength 
observed at the LHC at the expense of introducing a tension with EWPD and/or a low 
cutoff due to perturbativity and vacuum stability. In the models with extra vector-like 
leptons, Fermi gamma-ray line can give new bounds on the masses of the extra leptons, 
which constrain the region of the enhanced Higgs-to-diphoton rate. When Fermi gamma- 
ray line is explained by the two-photon line, which is the case in all the minimal models, 
the lighter charged lepton must be heavier than 130 GeV, which in turn bounds the ratio 
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of Higgs-to-diphoton rate to the SM value to at most 1.4(1.5) even for the low cutoff of 
10(1) TeV. 
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Appendix A: Axion couplings 

From eq. (J1J with eqs. f|T6|) and ffl~9]) . in the basis of mass eigenstates, we obtain the 
interaction terms for the lighter axion and the lighter CP-even Higgs in terms of Weyl 
spinors as 

- C int = i^-axx + i — a^4^4 + ^a(Ai/i/i + A2/2/2 + A3/1/2 + A4/2/1) 

v s v s 

-h {yifji + 2/2/2/2 + 2/3/1/2 + 2/4/2/1) + h.c. (A.l) 

where a is the mixing angle between two CP-even Higgs bosons and the Higgs couplings 
for all models that we considered in the text are 
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and the axion couplings are: for Model I 
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(A.6) 
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(A.7) 




(A.8) 




(A.9) 




(A.10) 



mi sin Of 

A 2 = (mi sinOf + m x cos 6 A, (A. 11) 

v s m h 

mi cos Or . 

A3 = -(mz sinflf + m.j; cos tM, (A. 12) 

sin Of . . , . 

A 4 = — (m^ cos 67 — m/ sin Of). (A. 13) 

v s m fl 

In the decoupling limit of the extra Higgs fields, we need to impose a = (3 — | so that we 
get — sin a/v d = cos a/v u = 1/v. 

Appendix B: Gauge couplings to extra leptons 

For a pair of left-handed lepton doublet (^4, eA) T and right-handed singlet e% the elec- 
troweak interactions are 



£ gauge = ^^t(V 47 M(l- 7 5)e + h.c.; 



2^2 
2 sin cos W 



1/4, e 

-e^e 7 M e (B.i; 



where e = (e 4 , e^) T , and 



1 , 

= a V4 = 2' ( B - 2 ) 

1 1 

v e = -- + 2 sin 2 6V, a e = --. (B.3) 

Thus, for a vector-like lepton doublet and a vector-like lepton singlet, we need to include the 
same electroweak interactions for the vector-like partners, a right-handed lepton doublet 
(£4,64) and a left-handed singlet e 4 . Then, we get the electroweak interactions for them 
in terms of Dirac spinors as 

£g„ug, = ^=^(i>T"(l-T 5 )e + P7"(l+7 5 )e' + hx.) 

e,e' 

-eA^erfe + e^e) (B.4) 

where v = (v^v^ T , e' = (el,e\) T , and the neutral current interactions are v e > = v e and 
a e = —a e i — — ~. Therefore, using the rotation matrices in eqs. (TTBT) and (IT5|) . the above 
interactions become in the basis of mass eigenstates, 
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2 sin 6w cos 6 W M L 
+ \a e {V^V u - V^V 2i )F^{l - 7 5 )F, + \a e (\j) x U u - U] 2 U 2l )F^(l + 7 5 )^ 
-A,(F 1 YF 1 + F 2 YF 2 ) (B.5) 

where Fi = f\) T and F 2 = (f 2 , f\) T are Dirac spinors with masses m/ i; m/ 2 . 

Appendix C: Matrix elements for the axion decay 

In this section, we present the expressions for the matrix elements for the axion decay. 
Suppose that the vector-like fermions, ipi = (f\,f 2 ,i>), have the following axion and 
gauge interactions, 

C f = iX^a 5 ^ + ^{gv t] - g'v^ 5 )^. (C.i) 

We consider the case that g' Vi . = and A^, gw i:j and g ltj are diagonal as for the simple 
charged fermion mass matrix with mi = m e and m x = m x in Model I. So, we denote the 
diagonal components by A^ = X^ u , gy i = gy tl and g 7i = g liV Then, the matrix elements 
for the decay channels of the axion are M. = e* fl (ki)e* u (k 2 )A4^ u , with 

M, v (a 77) = ^ E ^ Al N> ( C - 2 ) 

A ^ A ^ 

M^a^ZZ) = ^e, vp X 2 K\Y J J ^Mn,pf) + J ^A 2 (T V) p z u ) 



% 

MtuP i + ^ML MT2 , pl _ (c . 3) 

AV(« -> Z 7 ) = ^iW^E^^^^^f), (C.4) 

AU« WW) = W ^ [ E ^ ^ ^ , -)+E ¥^ Mn, P?, —) 
4:7T Z \-. 2m u m v ^ m fi m fc * 

(C.5) 



where ti j2 = 4m ^ Jm\, r v = Ami /ml, an d p\ :2 = 4m^ J my, pX = 4m^/— 2 



v • 



Ai(r) = 2/ da;/ — — - = r arcsin 2 (l/Vr), (C.6) 
Jo Jo A i( T ) 

A 2 (r,/) = 2 [ dx [ X dy A / v , , (C.7) 
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A 3 (r,p v ,b) 



1 

dx 
u Jo 
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A 4 (r,p 
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dx 
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(C.8) 
(C.9) 
(CIO) 

(C.ll) 
(C.12) 



4 

Ai(r) -^v(x + y)(l-x-y), 
Ai(r) - (x + y)(l - x - y) - (1 - 6 2 )(1 - x - y), (C.13) 
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Ai(r) - — ^ 2/(1 - x- y), 
4 

Ai(r) -x-y). 



P* 



(C.14) 
(C.15) 



Then, the axion decay rates are 
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with = gy./An, etc. When the axion is kinematically allowed to decay into a dark 
matter pair, that is, for m a > 2m x , we need to include the extra channel with decay rate 
given by 

|AJ 2 / 4m 2 \ 1/2 
r a (a^xx)= l -TF L rn a (l-^) . (C.20) 



167T 



Furthermore, the amplitudes for the DM annihilation into a pair of gauge bosons, V\ 
and V 2 , are 



M 



xx-m 



s-m 2 a - iT a m a 



v 2 ) , 



(C.21; 



where m 2 , T{ in M. xx ^ a and M.a-^v x v 2 are replaced by s and 4m 2 /s, respectively, in the 
axion decay amplitude part. As a consequence, the annihilation cross sections are given 
by 
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Appendix D: Renormalization group equations with 
extra leptons 



In the presence of the Yukawa coupling of extra vector-like leptons to the Higgs, we 
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consider the renormalization group (RG) equations below the masses of heavy Higgs states 
in the decoupling limit. 

Following the results [27], we obtain the RGEs, -4£ = /3 Pi , with t = \n(fi/m t ), with a 
vector-like doublet fermion and a vector-like singlet fermion, for the SM gauge couplings, 
the top Yukawa coupling, the extra Yukawa couplings and the Higgs quartic coupling, 

16tt 2 /V = 1V 3 , 16-% = -^ 3 , 16vr 2 /3 53 = -7gl (D.l) 



i6vr% 
167r% 
167r% 



f 2 i 2 , -2 -L'/2 ^2 o 

y* +Vi +Vi ~J^9 - 8 

I/J + Vi + 3 ^ ~ _ \9 

~ ( ~2 2 o 2 /2 ^ 

2/^2% +Vi +3 Vt ~-±9 



(D.2) 
(D.3) 
(D.4) 



16tt 2 /3 a = 24A 2 -(3^ 2 + V)A + ^(^ + 2^V + 3/) 

+4(3y 2 + y 2 + y 2 )A - 2(3 % 4 + yf + yf). (D.5) 

We note that we have ignored the quartic couplings of the singlet mediator 5 in the above 
RGEs. The Higgs-singlet coupling can contribute a positive term to the RG equation for 
the Higgs quartic coupling so it can help improve the vacuum stability a bit. 

On the other hand, for a vector-like doublet fermion and a triplet fermion, we get the 
corresponding RGEs, 

45 7 
16vr 2 /V = —g 3 , 16vr% = --g\ 16tt 2 /3 S3 = -7gl (D.6) 



16-% = yt^ + lyf + lyf-^-lg'-Sgl), (D.7) 

16-% = y^yf + yf + Sy^-^-^g 2 ), (D.8) 

16-% = m^yl + yf + ^-lgV-^g 2 ), (D.9) 

167r% = 24A 2 -(3^ 2 + V)A + ^(^ + 2^V + 3/) 

+2(6y 2 + 3yf + 3yf)\ - fa + \y\ + \yf) . (D.10) 
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